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We argue that the nonrelativistic Hamiltonian of p x + ip y superconductor in two dimensions can be 
derived from the relativistic Jackiw-Rossi model by taking the limit of large Zeeman magnetic field 
and chemical potential. In particular, the existence of a fermion zero mode bound to a vortex in the 
Px + ip y superconductor can be understood as a remnant of that in the Jackiw-Rossi model. In three 
dimensions, the nonrelativistic limit of the Jackiw-Rebbi model leads to a "p + is" superconductor 
in which spin-triplet p-wave and spin-singlet s-wave pairings coexist. The resulting Hamiltonian 
supports a fermion zero mode when the pairing gaps form a hedgehoglike structure. Our find- 
ings provide a unified view of fermion zero modes in relativistic (Dirac-type) and nonrelativistic 
(Schrodinger-type) superconductors. 
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I. INTRODUCTION 

Fermion zero modes bound to topological defects have 
been discovered by Jackiw and Rebbi in 1976 (Ref. Q]) 
and recently received renewed interest in condensed mat- 
ter physics (see, for example, Ref. y). Vortices in a cer- 
tain class of superconductors in two dimensions (2D) sup- 
port zero-energy Majorana bound states and obey non- 
Abelian statistics 3 , which can be potentially used for 
topological quantum computation^. Although vortices 
in the ordinary nonrelativistic s-wave superconductor do 
not support Majorana zero modes, the weakly paired 
phase of the p x + ip y superconductor, which is believed 
to be realized in S^RuO^, does support Majorana zero 
modes bound to vortex cores 3 ^^. 

R is also known from the pioneering work by Jackiw 
and Rossi that the relativistic s-wave superconductor in 
2D (Jackiw-Rossi model) has similar properties 8 . Re- 
markably it has been shown that such a system can be re- 
alized on the surface of the three-dimensional (3D) topo- 
logical insulator in contact with the s-wave superconduc- 
tor—. Besides these examples, there is a number of pro- 
posals to realize Majorana zero modes using heterostruc- 
tures of semiconductor and superconductor— ~—, super- 
conductor and ferromagnet 13 , and quantum (anomalous) 
Hall state and superconductor 14 . 

Although the nonrelativistic p x + ip y superconduc- 
tor and the relativistic Jackiw-Rossi model share similar 
properties, the existence of a fermion zero mode bound 
to a vortex has been discussed separately in the two sys- 
tems£&i£~— . In this paper (Sec. [n}, we argue that they 
are actually linked by showing that the former Hamilto- 
nian can be derived from the latter by taking the limit 
of large Zeeman magnetic field and chemical potential. 
In particular, the fermion zero mode bound to a vortex 
persists under taking this limit. 

Then in Sec. Mil we turn to the relativistic Jackiw- 



Rebbi model in 3D, which is known to exhibit a fermion 
zero mode associated with a pointlike topological defect 
(hedgehog) 1 ' 18 ' 20 ' 21 . The limit of large mass and chemical 
potential (nonrelativistic limit) leads to a "p-Ms" super- 
conductor in which spin-triplet p-wave and spin-singlet 
s-wave pairings coexist. We show that the resulting non- 
relativistic Hamiltonian supports a fermion zero mode 
when the pairing gaps form a hedgehoglike structure. 



We note that the analysis presented in this paper is 
largely motivated by the recent paper by Silaev and 
Volovik— : The nonrelativistic Hamiltonian of the Balian- 
Werthamer (BW) state of the superfluid 3 He was derived 
from the relativistic superconductor with the odd parity 
pairing^ 3 - and their topological properties were studied. 
In this paper, we shall broadly use "relativistic" to indi- 
cate Dirac-type Hamiltonians and "nonrelativistic" to in- 
dicate Schrodinger-type Hamiltonians. For readers' con- 
venience, references to the main results are summarized 
in Table |U 



TABLE I: References to the equations in which the Hamilto- 
nian, zero-energy solution bound to a defect, and its normal- 
izability condition are shown for the relativistic model and its 
nonrelativistic descendant both in 2D and 3D. 
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II. JACKIW-ROSSI MODEL IN 2D 
AND ITS NONRELATIVISTIC LIMIT 

A. Jackiw-Rossi model and fermion zero mode at a 
vortex 

We start with the Hamiltonian describing 2D Dirac 
fcrmions coupled with an s-wave pairing gap (Jackiw- 
Rossi 8 or Fu-Kane^ model) 



H =- I <lx¥m> 



(1) 



If we work in polar coordinates (r, 9) with the gap func- 
tion given by the vortex form 

A(x,y) = \A(r)\e m0 with |A(oo)| > 0, (5) 

it is easy to find the explicit zero-energy solution for odd 
N w = n (Ref. H 
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(2) 



This Hamiltonian can be realized on the surface of the 
3D topological insulator in contact with the s-wave su- 
perconductor 9 , h is the Zeeman magnetic field and p 
is the chemical potential. When the pairing gap A is 
spatially dependent, p = (p x ,Py) has to be regarded as 
derivative operators (— id x , —id y ). The energy eigenvalue 
problem is 



(3) 



When h and p are both zero, the number of fermion 
zero modes (e = 0) bound to a vortex formed by 
A(x,y) = A1-HA2 is determined by the winding number 
of the two scalar fields 8 ' 15,18 
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(4) 



where A a = A a /yAf + Af and the line integral is taken 
at spatial infinity. However, in the presence of h and /i, 
the index theorem is no longer valid: h and \i terms in 
the Hamiltonian can couple zero modes and they become 
nonzero energy states so that two states form a pair with 
opposite energies. Therefore, in general, only one zero 
mode survives for odd N w while no zero mode survives 
for even Nj^-. a 



Two exceptional cases are h = ±/i. When N w > (< 0), the 
zero-energy solutions at h = fi = still solve Eq. l|3} with 
h = +(■—)/! ^ and thus there are \N W \ zero modes. This 
can be easily seen if one rewrites the Hamiltonian in the ba- 
sis where the chiral operator defined in Eq. (12 1 f) has the form 
X = diag(l, 1, — f, — f) and recognizes that the \N W \ zero-energy 
solutions are eigenstates of \ with the eigenvalue +1 (—1) for 
N w >0(<0). 



and v\ = — u*, v 2 — u\ with an integer I = (n— l)/2. We 
note that the zero-energy solution, Eq. (J6)), is normaliz- 
able as long as 



^ 2 + |A(oo)| 2 >h 2 



(7) 



is satisfied and there is a topological phase transition at 
l? + |A| 2 = h 2 [see also Eq. flU} below]. 



B. Derivation of p x + ip y superconductor and 
fermion zero mode 



We now derive a clear connection between the Jackiw- 
Rossi model and the nonrelativistic p x +ipy superconduc- 
tor. Suppose we are interested in the low-energy spec- 
trum of Hamiltonian ((5J in the limit where both h > 
and /i > are equally large 



|vV + |A| 2 -/i|«/i~m- 



(8) 



The low-energy spectrum in such a limit can be obtained 
by eliminating small components u 2 and v<j&. Substitut- 
ing the following two equations from Eq. ([3j): 



(e + h + /i) u 2 = p + ui + Avi 
(s — h - fi) v 2 — —p-Vi + A*m 

into the remaining two equations, we obtain 

p 2 ui+ p_Avi —Ap_vi + \A\ 2 ui 

e + h + fi e — h — p 

p 2 vi - p+A* Ul A* P+Ul + \A\ 2 Vl 



(9) 



{e — h 4- p) u\ 
(e + h — p) vi 



h — p, 



/' 



(10) 



Here we introduced p± = p x ± ip y . 

In the limit under consideration, Eq. (J8J), we can ne- 



glect £ compared to h + p and approximate v/ ' p 2 + \A\ 2 
by h. The remaining components u\ and v\ obey the new 
energy eigenvalue problem 



Mnr i{i3_,A nr } 

2 



V 

2in 




(11) 



where we defined the nonrelativistic mass, chemical po- 
tential, and pairing gap as 



m = h, ^„ r = v V + |A| 2 - h 
The resulting Hamiltonian 
( *- 



A nr = -. (12) 
h 



/i nr ±{p_,A nr } 



(13) 



describes the nonrelativistic p x + ip y superconductor. We 
note that when h < 0, one obtains the Hamiltonian of the 
Px — T-Py superconductor where p+ and p- are exchanged 
in Eq. (H3J). 

The first nontrivial check of this correspondence is the 
comparison of spectrum in a uniform space where A is 
constant. The relativistic Hamiltonian (2) has the energy 
eigenvalues 



e 2 =p 2 + h 2 + p 2 - 



|A| 2 ±2v/pV + ^ 2 (M 2 + |A| 2 ). (14) 



Its low-energy branch (lower sign) at small p is correctly 
reproduced by the energy eigenvalue of the nonrelativistic 
Hamiltonian CO 



P 
2m 



V|A B 



(15) 



under the assumptions in Eq. (j8|). 

Because the above "nonrelativistic limit" does not rely 
on the spatial independence of A, the fermion zero mode 
found in Eq. (5]) persists into the p x +ip y superconductor, 
Eq. (|13p . In order to demonstrate it, we consider the 
simplified vortex configuration with a constant |A nr | > 

A nr (x,y) = e me \A nr \. (16) 

When n is odd, we can find the explicit zero-energy so- 
lution (e = 0) to Eq. (HU) (Ref. [13) 



Ml = J t 



2m / u„ r - (m|A nr |) r 



■^i+iW— m\ A nr |r /ir-^ 



and V\ = —u\. One can see that this zero-energy solu- 
tion is the direct consequence of that in Eq. ([6]) because 
Eqs. © and (T2]) lead to 



|A„ 



(i —h=(n+h){n—h)tt 2m ( fi 

'(18) 

Thus we have established that the existence of a 
fermion zero mode bound to a vortex in the p x + ip y 
superconductor, Eq. (|13p . is a remnant of that in the 
Jackiw- Rossi model, Eq. (2). In particular, the condi- 
tion for the normalizability of the zero-energy solution, 
Eq. ([7]), is translated into 



/inr > 



(19) 



TABLE II: Properties under the time-reversal operator T. r- 
matrices act on the particle-hole space and o (x) indicates 
even (odd) under T. Replacement of to by tz exchanges the 
roles of Ai and Ag. 



T 



-T/T 



CT2 <X>T 
01 ®Tl 



-1 

+ 1 
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which coincides with the well-known topological phase 
transition in the p x +ip y superconductor existing at /i nr = 
(Refs. drill) [see also Eq. (25) below]. Our finding 
also clarifies why a vortex with winding number N w in 
the p x + ip y superconductor cannot support |iV,„| zero 
modes in contrast to in the Jackiw-Rossi model with h = 
H = 16 i 17 . In order to derive the p x + ip y superconductor 
as a nonrelativistic limit of the Jackiw-Rossi model, one 
needs to introduce h and /z which split an even number 
of zero modes into positive- and negative-energy states. 
Therefore, only one zero mode survives for odd N w in the 
P x + ipy superconductor. 



C. Altland-Zirnbauer symmetry class (Refs. [27|j28h 
and topological invariant 

Finally, we note the Altland-Zirnbauer symmetry class 
of the Hamiltonians that we have investigated in this sec- 
tion. Hamiltonian (2J) with spatially dependent A x 2 7^ 
has the charge conjugation symmetry 



C^HC = -H* 



with 



C = 




(20) 



The properties of each term under the time-reversal op- 
erator T (T~ 1 'HT — H* at h = fi = Ai, 2 = 0) are sum- 
marized in Table [Til In particular, the so-called chiral 
symmetry, 



X - 1 H X = -H with x = r 3 ^ 



(21) 



is present only if h = fi = and essential for the index 
theorem, Eq. (4]). Therefore, the relativistic Hamiltonian 
(2|) with /i,/i^0 and thus the resulting nonrelativistic 
Hamiltonian (TT31) belong to the symmetry class D. 

According to Refs.l23andl30l. the class D Hamiltonians 
defined in compact 2D momentum spaces can be classi- 
fied by an integer- valued topological invariant, which is 
the first Chcrn number 31-33 



with 



r - — Id ( 9a,y - ® ax 

1 ~ 2-k J \dp x dp y 



o»(p) = Yl *! £o ' p l a~ l £a ' p " 



£a<0 



dpi 



(22) 



(23) 



We shall use Eqs. (|22j) and (123)) as a definition of the 
topological invariant C\ even for relativistic (Dirac-type) 
Hamiltonians while C\ in this case can be a half integer. 
However, for superconductors, C\ is always an integer 
because of the Nambu-Gor'kov doubling. We find that 
the topological invariant for the relativistic Hamiltonian 
^ is given by 



Ch 




-sgn(/i) 



for M 2 + |A| 2 > h 2 



for 



/<-" 



\A\ 2 <h 2 , 



(24) 



while the topological invariant for the nonrelativistic 
Hamiltonian (| 13[) is given by 



Ci = 



for p, nY > 
for /j nr < 0. 



(25) 



Therefore, in general, the topological invariant of the mo- 
mentum space Hamiltonian is not preserved by the non- 
relativistic limit. b 

Nevertheless, both values of C\ computed for the rel- 
ativistic and nonrelativistic Hamiltonians are consistent 
with recent conjectures relating the topological invari- 
ant of a momentum space Hamiltonian to the number of 
fermion zero modes bound to a vortex 2 ^. For class D 
superconductors defined in compact momentum spaces 
(as is the case for nonrelativistic Hamiltonians), Teo and 
Kane in Ref. l2| conjecture that the number of fermion 
zero modes is 



CiN w mod 2. 



(26) 



This formula gives v = 1 for /i nr > and v — for 
/i nr < for an odd winding number N w . On the other 
hand, Santos et al. in Ref. 34 do not constrain Hamiltoni- 
ans to be defined in compact momentum spaces, allowing 
for relativistic (Dirac-type) Hamiltonians, and conjecture 
that the number of fermion zero modes is 



v = (d+N f )N w mod 2, 



(27) 



where Nf is the number of Dirac flavors [Nf — 1 for the 
Jackiw- Rossi model, Eq. @, and Nf ~ for the p x + ip y 
superconductor, Eq. (|13[l ]. For an odd winding number 
N w , their formula gives v = 1 for p 2 + |A| 2 > h 2 and 
fj, nT > and v = for /j 2 + |A| 2 < h 2 and p, nr < 0. There- 
fore, the conjectured counting of fermion zero modes in 
terms of the momentum space topological invariant works 
both in the relativistic and nonrelativistic Hamiltonians, 
even though the value of C\ is not preserved by the non- 
relativistic limit. 



b The topological invariant can be matched if we properly regu- 
larize the large p behavior of the relativistic Hamiltonian: Re- 
placing h in Eq. ([2j by h (l + |e|p 2 ) , the Chern number becomes 
sgn(h) for fi 2 + |A| 2 > h 2 and for fi 2 + |A| 2 < h 2 , which 
coincides with that of Eq. 113t . 



III. JACKIW-REBBI MODEL IN 3D 
AND ITS NONRELATIVISTIC LIMIT 

A. Jackiw-Rebbi model and fermion zero mode at 
a hedgehog 

In this section, we extend the above developed anal- 
ysis to three dimensions. For this purpose, we consider 
the following Hamiltonian describing 3D Dirac fermions 
coupled with three real scalar fields (A = Ai + iA2 and 
A 8 ): 



H 



1 



dx &W& 



(28) 



with vE'' = (ip\ —iip T a2) and 



H 



a ■ p + (3m — [i — «7 5 / 9A3 

A* 



a. ■ p + /3m + jj, + ?7 5 /3A3 
(29) 



This Hamiltonian with zero mass m and zero chemical 
potential fj,, after an appropriate unitary transformation 
and renamings {(3 f-> «7 5 /?, Ai — > <j)\, A 2 — > — fa, A 3 — >• 
—fa), was studied initially by Jackiw and Rebbi^ and re- 
cently by Teo and Kane 21 in the context of ordinary and 
topological insulators coexisting with superconductivity. 
When the scalar fields Ai i2 ,3 are spatially dependent, 
p = {p x ,Py,p z ) has to be regarded as derivative opera- 



tors (- 

is 



id x , —id y , —id z ). The energy eigenvalue problem 



(30) 



where ui,2 and v\^ are two-component fields. Here we 
employ the standard representation of Dirac matrices 
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U/ 



a- " 
\a 0/ 

and hence 
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z 7 5 /3 = 



-it 
M , 



(31) 



(32) 



When to and /i are both zero, the number of fermion 
zero modes (e = 0) bound to a hedgehog formed by 
Ai,2, 3(^,2/, z) is determined by the winding number of 
the three scalar fields ^ 18 ' 20 

Index'H = — l dS t e ijk e a bc^adjA b d k A c = N w , (33) 



where A a = A a /-\/A 2 + A2 + A3 and the surface inte- 
gral is taken at spatial infinity. However, in the presence 



of m and fi, the index theorem is no longer valid: m and fi 
terms in the Hamiltonian can couple zero modes and they 
become nonzero energy states so that two states form a 
pair with opposite energies. Therefore, in general, only 
one zero mode survives for odd N w while no zero mode 
survives for even N^r^. 

Here, instead of the symmetric hedgehog (Aj oc ajj), 
we assume the hedgehoglike configuration in which A 1.2 
depend only on (x, y) and form a vortex and A3 depends 
only on z and forms a kink. They have the same winding 
number but the latter has the advantage that an analytic 
solution can be found even with m, /j ^ 0. If we work 
in cylindrical coordinates (r, 8, z) with the gap functions 
given by the forms 



and 



A{x,y) = \A(r)\e in6 with |A(oo)| > (34) 



A 3 (z^±w)^±|A 3 |, (35) 



it is easy to find the explicit zero-energy solution for odd 
N w = n (Ref. \M) 



y/ \i + m Ji f v /i 2 — m 2 r)e 4 1 



y/fj, — mJi + \ ( yj [i 2 — m 2 r) e^ t+t 
x e iW-f r dr'\A(r')\-f dz'A 3 (z') 



(36) 



and v\ = i<j2Ui, v^ — 102112 with an integer I = (n— l)/2. 
On the other hand, when 



A 3 (z^±oo)^t|A 3 



(37) 



with the same A(x, y) in Eq. (I34[) . we have N w = —n and 
the zero-energy solution in Eq. (J36]) is replaced by 



\i+i8 



x e 



yjii + m Ji+i (y/fi 2 - m 2 r) 
y/jjL — m Ji I \/fj, 2 — m 2 r) e~ 5 ' 


iW~f r dr'\A(r')\+J z dz'A 3 (z') 



(38) 



We note that the zero-energy solution, Eq. (|36l) or (13"8"1) . 
is normalizable as long as 



/1 2 + |A(oo)| 2 > m 2 



(39) 



is satisfied. 



Two exceptional cases are m = ±/7. When N w > (< 0), the 
zero-energy solutions at m = \x = still solve Eq. I|30| l with 
m = +(— )fJ. 7^ and thus there are \N W \ zero modes. This 
can be easily seen if one rewrites the Hamiltonian in the ba- 
sis where the chiral operator defined in Eq. {55} has the form 
X = diag(l, 1, — 1, — 1) and recognizes that the \N W \ zero- 
energy solutions are eigenstates of \ with the eigenvalue +1 (—1) 
for N w > 0(<0). 



B. Derivation of p + is superconductor and fermion 
zero mode 



We now study the nonrelativistic limit of the above 
Jackiw-Rebbi model. Suppose we are interested in the 
low-energy spectrum of Hamiltonian (j2"9")l in the limit 
where both m > and fi > are equally large 



e, \y/n 2 + \A\ 2 -m\ < 



777 ~ jj,. 



(40) 



The low-energy spectrum in such a limit can be obtained 
by eliminating small components u 2 and v?&- Substitut- 
ing the following two equations from Eq. ([30| : 



(41) 



(e + m + jx) u 2 = (a- ■ p - iA 3 ) u\ + Av\ 
(e - m — fi) v 2 = - (a- ■ p + iA 3 ) V\ + A*ui 

into the remaining two equations, we obtain 

(e — 777 + /j) Ml 

_ [p 2 + A§ - o- • (9A 3 )] ui + (cr -p + iA 3 ) Atji 
e + m + fj, 

-A(cr-p + 7:A 3 )7;i + |A| 2 ui 



e — m — [i 
(e + 777 — li) Vi 

_ [p 2 + A 2 + a • (0A 3 )] v x - {a ■ p - iA 3 ) A* Ul 



£ — 777 — \X 

A*(CT-P-7A 3 )771 + |A| 2 7; 1 
£ + 777 + Li 



(42) 



Here the derivative operator d in er • (<?A 3 ) acts only on 
A3. 

In the limit under consideration, Eq. (|40p . we can ne- 
glect £ compared to m + /1 and approximate y/fi 2 + |A| 2 
by 777. The remaining components ui and v\ obey the 
new energy eigenvalue problem 



Ml 



~^§^ \{a-p,At} + iK 



p 

2 X" Pi t±t I llS s 2m + ^ nr 2m 



(43) 



where we defined the nonrelativistic chemical potential 

as 



/x nI = x/J^+jAj 



A 2 
2m 



(44) 



and the spin-triplet p-wave and spin-singlet s-wave pair- 
ing gaps as 



A , A A 3 A 

A t = — and A s = . 

777 777 



(45) 



The resulting Hamiltonian 



rl m — 



p 

2m 



g-(8A 3 



h{<r-P,*t} 



2m 

-iA* 



\{a ■ p, A t } + iA s 



2m ' A*n 



■(flA 3 ) 



2m 



describes the p + is superconductor in which spin-triplet 
p-wave and spin-singlet s-wave pairings coexist. A s can 
be complex but its phase is locked to the phase of A t 
[see Eq. (|15j)] and thus there are three independent de- 
grees of freedom. The last term in the diagonal elements 
resembles the Zeeman coupling er ■ B with "magnetic 
field" Bi = —diA^/(2ra) generated by the gradient of 
A3 = A s / At. We note that the nonrelativistic Hamil- 
tonian Pr?|) in the absence of A s is the BW state of the 
superfluid 3 He and studied in Ref. |2J- 

The first nontrivial check of this correspondence is the 
comparison of spectrum in a uniform space where A and 
A3 are constant. The relativistic Hamiltonian (J2T))) has 
the energy eigenvalues 



e 2 = p 2 



/'- 



A 2 



±2Jp 2 p 2 +m 2 (fi 2 + \A\ 



M 2 A 2 . 



(47) 



Its low-energy branch (lower sign) at small p and A3 
is correctly reproduced by the energy eigenvalue of the 
nonrelativistic Hamiltonian (1461) 



P 

1m 



fh: 



V|A t | 



|A S 



(48) 



under the assumptions in Eq. (|4U|) . 

Because the above nonrelativistic limit does not rely 
on the spatial independence of A and A3, the fermion 
zero mode found in Eq. (l36l) or (|38|) persists into the 
p + is superconductor, Eq. (|4"rj)) . In order to demonstrate 
it, we consider the simplified hedgehoglike configuration 
resulting from Eqs. (l34l) . (l35l) . and (j37|) with constant 
I A t I > and |A S | > 0: 



A t (x,y) = e me \A t \ and 
A s (x,y,z)=±e m9 sgn(z)\A s 



(49) 



When n is odd, we can find the explicit zero-energy so- 
lution (e = 0) to Eq. 03) 



»i 



Ji 



2mp m — (m\At 







(50) 



corresponding to the upper sign in Eq. (|49[) , or 


J1+1 



«1 



'2rof* nr -(ro|At|)* + 

x e ft+*(J+l)9-Tn|A t |r-|^.||z| 



(51) 



TABLE III: Properties under the time-reversal operator T. 
r-matrices act on the particle-hole space and o (x) indicates 
even (odd) under T . Replacement of to by T3 exchanges the 
roles of Ai and A2. 
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7^/r 
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A3 


Ai 


A 2 
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-1 





X 
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7 5 Q 2 ® To 


-1 








X 
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/3 OL2 ® T\ 


+1 
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X 











7 5 /? q 2 ® n 


-1 


X 





X 









corresponding to the lower sign in Eq. (1491) . and v\ = 
ia2u\. One can see that this zero-energy solution is the 
direct consequence of that in Eq. (IM)) or (1551) because 
Eqs. PD ]) , (jit ]) , aad (|55 ]l lead to 



2 2 
p — m 



2to 



/'-!) 



: 2m\A t 



(52) 



Thus we have established that the existence of a 
fermion zero mode bound to the hedgehoglike structure, 
Eq. (|4^|) . formed by A t and A s /A t in the p + is super- 
conductor, Eq. (I46p . is a remnant of that in the Jackiw- 
Rebbi model, Eq. (|2"9")h In particular, the condition for 
the normalizability of the zero-energy solution, Eq. (|39p , 
is translated into 



A*nr + 



2m 



A, 



>0. 



(53) 



C. Altland-Zirnbauer symmetry class (Refs. I27lj28h 

Finally, we note the Altland-Zirnbauer symmetry class 
of the Hamiltonians that we have investigated in this sec- 
tion. Hamiltonian (|2"9")) with spatially dependent Ai^.3 ^ 
has the charge conjugation symmetry 



c^nc = -n* 



with 



C = 




(54) 



The properties of each term under the time-reversal op- 
erator T {T^HT = rl* at m = p = Ai, 2 ,3 = 0) are 
summarized in Table Hill In particular, the so-called chi- 
ral symmetry, 



x - 1 n x = -n 



with 



X = 




(55) 



is present only if m = \x = and essential for the index 
theorem, Eq. (|3"3")1 . Therefore, the relativistic Hamilto- 
nian (|29[) with m, p ^ and thus the resulting nonrel- 
ativistic Hamiltonian (1461) belong to the symmetry class 
D. There is no topological classification of class D Hamil- 
tonians in 3D momentum spaces 29,30 . 



IV. SUMMARY 

Wc have studied the nonrelativistic limit of the Jackiw- 
Rossi model in 2D and the Jackiw-Rebbi model in 3D, 
both of which are known to exhibit fermion zero modes 
associated with pointlike topological defects (vortex and 
hedgehog). We showed that the nonrelativistic limit of 
the 2D Jackiw-Rossi model leads to the p x + ip y super- 
conductor. Because the fermion zero mode persists under 
taking this limit, we obtain a clear understanding of the 
existence of a fermion zero mode bound to a vortex in 
the p x + ip y superconductor as a remnant of that in the 
Jackiw-Rossi model. Similarly, the nonrelativistic limit 
of the 3D Jackiw-Rebbi model leads to the p + is super- 
conductor in which the spin-triplet p-wave pairing gap 
At and the spin-singlet s-wave pairing gap A s coexist. 
We showed that the resulting Hamiltonian supports a 
fermion zero mode when A t and A s /A t form a hedge- 
hoglike structure. Fermion zero modes in the supercon- 
ductors studied in this paper correspond to Majorana 



fermions and the associated pointlike defects obey non- 
Abelian statistics both in 2D (Refs. [30) and 3D2L2S. 

Our findings provide a unified view of Majorana zero 
modes in relativistic (Dirac-type) and nonrelativistic 
(Schrodinger-type) superconductors. It should be pos- 
sible to generalize our analysis to other interesting cases 
and find new examples of nonrelativistic Hamiltonians, 
which are more common in condensed matter systems, 
with topological properties that descend from Dirac-type 
Hamiltonians, which are generally easier to analyze. 
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